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We compute the light and strange quark masses me = (m u + m<j)/2 and m 3 , respectively, in 
full lattice QCD with Nf — 2 flavors of light dynamical quarks. The renormalization constants, 
which convert bare quark masses into renormalized quark masses, are computed nonperturbatively, 
including the effect of quark-line disconnected diagrams. We obtain mf s (2GeV) = 4.7(2) (3) MeV 
and mf s (2GeV) = 119(5)(8) MeV. 
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PACS numbers: 12.15.Ff,12.38.Gc,14.65.Bt 
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The light and strange quark masses are among the least well known parameters of the Standard Model. The reason 
is that quarks are confined, so that the masses must be determined indirectly through their influence on hadronic 
observables. This requires nonperturbative techniques. One such technique is lattice QCD. 

The quark masses obtained directly in lattice calculations are bare quark masses at the cut-off scale a~ 4 , where a 
Qh, denotes the lattice spacing. For the lattice numbers to be useful for phenomenology, it is necessary to convert the bare 
quark masses to renormalized masses in some standard renormalization scheme. Because lattice perturbation theory 
converges badly, and the expansion coefficients are generally known to one loop order only, this ought to be done 
^-h nonperturbatively. In full QCD a one-loop perturbative renormalization of the mass operator is totally inadequate 
, even, as it does not account for the disconnected (flavor singlet) contribution shown in Fig. 1, which turns out to be 
comparable with the connected contribution at present lattice spacings. 

In quenched QCD, in which the effect of sea quarks is neglected (and hence quark-line disconnected fermion loops 
are absent), several groups |lj have carried out an entirely nonperturbative calculation of the light and strange quark 
masses. Remarkably consistent results have been found. Previous calculations in full QCD, both with Nf — 2 2] 
and Nf = 3 |3j flavors of sea quarks, employ perturbative renormalization to compute the relation between the bare 
and renormalized quark masses. These authors, except perhaps Eicker et al., find rather small values for the strange 
quark mass, which lie substantially below the central value quoted by the Particle Data Group [4|. 
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FIG. 1: Quark diagrams contributing to the renormalization of the mass operator (x). The left figure shows the connected 
(nonsinglet) contribution, the right figure the disconnected (singlet minus nonsinglet) contribution. Gluon lines have been 
omitted. 
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TABLE I: Hopping parameters of sea and valence quarks used in this calculation, together with their critical values. 
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5.20 


0.13420 
0.13500 
0.13550 


0.13340 0.13380 0.13420 0.13470 0.13530 0.13560 0.1360 0.13620 
0.13320 0.13370 0.13420 0.13450 0.13500 0.13530 0.1355 0.13570 
0.13320 0.13360 0.13400 0.13430 0.13480 0.13500 0.13530 0.13550 0.13570 


0.137550(49) 
0.136889(32) 
0.136457(23) 


0.136008(15) 


5.25 


0.13460 
0.13520 
0.13575 


0.13370 0.13400 0.13460 0.13490 0.13530 0.13550 0.13590 0.13610 
0.13370 0.13410 0.13450 0.13480 0.13520 0.13550 0.13580 0.13590 
0.13360 0.13390 0.13430 0.13460 0.13500 0.13520 0.13550 0.13575 0.13600 


0.137237(19) 
0.136883(13) 
0.136553(9) 


0.136250(7) 


5.29 


0.13400 
0.13500 
0.13550 


0.13400 0.13440 0.13490 0.13520 0.13550 0.13570 0.13590 0.13620 
0.13400 0.13430 0.13470 0.13500 0.13550 0.13570 0.13600 0.13610 
0.13390 0.13430 0.13460 0.13490 0.13530 0.13550 0.13580 0.13600 0.13630 


0.137516(33) 
0.137045(16) 
0.136816(11) 


0.136410(9) 


5.40 


0.13500 
0.13560 
0.13610 


0.13420 0.13480 0.13500 0.13530 0.13560 0.13590 0.13610 0.13630 0.13660 
0.13460 0.13500 0.13530 0.13560 0.13570 0.13595 0.13610 0.13650 
0.13470 0.13500 0.13530 0.13560 0.13580 0.13595 0.13610 0.13645 


0.137131(14) 
0.136966(12) 
0.136836(14) 


0.136690(22) 



In this Letter we shall present a first fully nonperturbative calculation of the light and strange quark masses in 
full QCD, including the effect of flavor singlet renormalization factors. We consider nonperturbatively 0(a) improved 
Wilson fermions with Nj = 2 flavors of degenerate dynamical quarks and the Wilson gauge field action. The calculation 
is done in two steps. We simulate dynamical gauge field configurations at four different values of the coupling, /3, 
and at three different sea quark masses each. The latter are specified by the hopping parameter K sca - The actual 
parameters, as well as the corresponding lattice spacings and pseudoscalar mass values, are shown in Fig. 2. We use 
the force parameter r$ to set the scale. On these configurations we then perform a partially quenched calculation of 
the pseudoscalar mass, allowing for different sea and valence quark masses, from which we derive the physical quark 
masses. In Table 1 we list the hopping parameters of valence (/-c V ai) and sea quarks considered in this calculation. 

The bare sea and valence quark masses are given by a,TO sca = l/(2/c soa ) — l/(2/Cg Ca ) and am va i = l/(2ft va i) — l/(2«g Ga ), 
respectively. We consider the case of degenerate valence quarks only. In Fig. 3 we show the partially quenched 
pseudoscalar mass mps(K sea , K va i)- The critical hopping parameter Kg ea is found by keeping (3 fixed and varying K sea 
until mps(K S ea,Ksca) = 0. Similarly, we introduce a critical hopping parameter of the valence quarks, «^ al , which 
is found by varying K va i until mps (ksch, «vai) — 0, while keeping f3, K sca fixed. Our calculation requires a precise 
determination of Kg Ca and k^bA- We perform a global fit of the form 

2 f 1 1 \ f 1 1 \ ( 1 1 
ar m PS = u [ - ) + V I ■ — I +w I — 

\ K val K sca / \ K sca K sca / -| 

\ ^sca ^sca / \ ^val ^sca, / 

to all 100 data points, where the parameters are taken to be second order polynomials in (3. The fit gave x 2 /dof = 0.67. 
The resulting values of «g ea and re£ al are given in Table 1. The force parameter ro/a was computed from the static 
potential. 

Taking the derivatives of the quark propagator with respect to am sea and am V ai, we obtain for the renormalized 
sea and valence quark masses 

m soa = Z^ n m sca , (2) 
m vai = 2m S (m val - m SGa ) + Z^m soa , (3) 

where and Z^ s are singlet and nonsinglet renormalization constants of the mass operator. Partially quenched 
chiral perturbation theory to NLO predicts 
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2 PS = [A + (B + C In mf al ) m*J m? al + (D + E In m« al ) (mf al ) 2 . (4) 



From |0J follows that m^ al vanishes where the partially quenched pseudoscalar mass vanishes, which happens at the 
value K V ai = K vai- ^ we insert this value into J2Jl we obtain the ratio 
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FIG. 2: Parameters of our dynamical gauge field configurations, together with lines of constant ro/a (solid lines) and lines of 
constant mpsro (dashed lines). The simulations are done on 24 3 48 (x) and 16 3 32 (O) lattices, respectively. 
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FIG. 3: The partially quenched pseudoscalar mass as a function of m va i at f3 = 5.40 for K sca = 0.1350 (O), 0.1356 (□) and 
0.1361 (A), together with the fit Q. The solid line and symbols refer to the case K va i = n sca . 
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FIG. 4: The ratio Z^/Z^ s at /3 = 5.20, 5.25, 5.29 and 5.40 (from top to bottom), together with a linear extrapolation to the 
chiral limit. 



0.6 



0.5 - 



a 

ISI 



0.4 - 



0.3 



1 1 1 

- CD 


1 1 1 


i i i i i 


i 1 1 1 


1 1 1 1 


- © 










" © 










- @ € 




eee o e 


e e e 


0" 


i i i 


i , , 


= 5.40, 

1 i i i i 


k« =0. 

sea 
1 i i i 


1356 " 

i i i i 







20 40 60 80 100 

(r fi) 2 



FIG. 5: The nonsinglet renormalization group invariant Z§ GI at (3 = 5.40, K sca = 0.1356 as a function of the renormalization 
scale fi, together with the fit to the plateau. 
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FIG. 6: The partially quenched pseudoscalar mass mps as a function of m^ al at /? = 5.40 for our three sea quark masses. The 
solid line shows the result of the fit for mfL = 0. 
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FIG. 7: The lattice spacing a extrapolated to the chiral limit. 
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In Fig. 4 we plot Z^/Z^ s for all data sets. The effect of the quark-line disconnected diagram is found to be significant. 
The numbers depend only mildly on the sea quark mass. 

It remains to determine Z^ s . We compute Zg S = (Z^ 5 )" 1 nonperturbatively in the RI- MOM scheme. The 
result is converted to the more popular MS and RGI schemes by a three-loop perturbative calculation ||. In Fig. 5 
we show the nonsinglet Zg GI as a function of the renormalization scale fi. We find that the nonperturbative scale 
dependence of Zg is matched by the three-loop conversion factor for (r fi) 2 > 20. We obtain Z§ GI from 

a fit to the plateau as indicated by the solid line. The result varies by a few percent only over our range of sea 
quark masses at any given f3 value. In the MS scheme at \i = 2 GeV we have [U Zf s (2GeV) = 1.461 Z§ GI . At 
our smallest lattice spacing, a k 0.07 fm, Zjf s « 0.6, which is certainly beyond the range of one-loop perturbation 
theory, tadpole- improved or not. 

Having unscrambled renormalized valence and sea quark masses, we are now able to fit our data by the partially 
quenched chiral formula and determine the physical quark masses from it. In the process we replace all masses 
m by dimensionless quantities mr^. It turns out that the data are not sensitive to logarithmic variations in the 
renormalized valence quark mass (parameters C and E). We therefore have chosen to fit our partially quenched 
pseudoscalar masses by the formula 

(m PS r f = [A + B mf ea r ] m^r + D (m^r ) 2 . (6) 

In Fig. 6 we plot our data for our largest (3 value. The slope of the data depends only rather weakly on the renormalized 
sea quark mass. Perhaps most of the effect is washed out by having used ro to set the scale. The solid curve shows 
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FIG. 8: The light and strange quark masses, together with the extrapolation to the continuum limit. The errors shown are 
statistical only. 
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the result of the fit in the limit of vanishing sea quark mass. We find good scaling properties. The fit parameter A 
varies by less than 5% over our range of /? values. 

To fix the scale r in physical units, we extrapolate recent dimcnsionless nucleon masses, m^r , found by the 
CP-PACS, JLQCD and QCDSF-UKQCD collaborations jointly to the physical pion mass, following [9j. This gives 
the value rg — 0.467 fm, which we will use here. A similar result was quoted in [10|. The average mass of the up and 
down quarks, mi = (m u + md)/2, is found from extrapolating mps to the physical ir° mass, setting m^ al = mf ca in (0. 
We obtain mf w {2 GeV) r = 0.00981(19), 0.00987(17), 0.00986(18) and 0.01044(19) at /3 = 5.20, 5.25, 5.29 and 5.40, 
respectively. Similarly, the strange quark mass, m s , is obtained from the lattice value of m^ al that brings mps to the 
physical K° mass, while m^ a is kept fixed at the corresponding physical sea quark mass mf. Owing to the fact that 
the valence quarks are degenerate, we then have mf = 2mf al — mf. This finally gives mf IS (2 GeV) r = 0.2525(50), 
0.2544(45), 0.2545(47) and 0.2671(48) at (3 = 5.20, 5.25, 5.29 and 5.40, respectively. 

To be able to extrapolate our results to the continuum limit, we need to know a/ro in the chiral limit. In Fig. 7 
we show our data for ro/a together with a fit (fitting function: exponential of a polynomial in (3, m sea ), and in Fig. 8 
we show the light and strange quark masses as a function of the chirally extrapolated lattice spacing. Because our 
fermionic action is nonperturbatively O(a) improved, we expect the error due to the finite cut-off to be at most of 
0(a 2 ). A linear extrapolation in (a/ro) 2 to the continuum limit is therefore appropriate. We estimate the systematic 
error on rg to be of the order of 7%. We then obtain 



where the first error is statistical and the second systematic. In particular, m s /mg — 26(1), in good agreement with 
leading order chiral perturbation theory. 

To summarize, we have performed a lattice calculation of the light and strange quark masses in full QCD with 
Nj — 2 flavors of light dynamical quarks. Our calculation differs from previous calculations in several respects. 
We use nonperturbatively 0(a) improved Wilson fermions and perform simulations at four different couplings with 
0.07 < a < 0.12 fm, which allows an extrapolation to the continuum limit. Furthermore, an entirely nonperturbative 
scheme of mass renormalization, both for sea and valence quark masses, is devised, including the effect of quark-line 
disconnected contributions. The identification of renormalized sea and valence quark masses greatly facilitates the 
extrapolation to the chiral limit. We believe that our method will be useful in other applications of partially quenched 
chiral perturbation theory as well. 

This work is supported by DFG under contract FOR 465 (Forschergruppe Gitter-Hadronen-Phanomenologie). The 
numerical calculations have been performed at LRZ (Miinchen) ||, EPCC (Edinburgh) DESY (Zeuthen) and 
NIC (Julich). 



[1] S. Aoki et ai, Phys. Rev. Lett. 82, 4392 (1999); J. Garden et al, Nucl. Phys. B 571, 237 (2000); M. Gockeler et ai, Phys. 

Rev. D 62, 054504 (2000); D. Becirevic, V. Lubicz and C. Tarantino, Phys. Lett. B 558, 69 (2003). 
[2] N. Eicker et ai, Phys. Lett. B 407, 290 (1997); D. Pleiter, Nucl. Phys. Proc. Suppl. 94, 265 (2001); A. Ali Khan et ai, 

Phys. Rev. Lett. 85, 4674 (2000) [Erratum-ibid. 90, 029902 (2003)]. 
[3] T. Kaneko et ai, Nucl. Phys. Proc. Suppl. 129, 188 (2004); C. Aubin et al, Phys. Rev. D 70, 031504 (2004). 
[4] S. Eidelma n et al. [ Particle Data Group Collaboration], Phys. Lett. B 592, 1 (2004); see also: R. Gupta, 

ai'Xiv:hcp-ph/0311033 
[5] M. Gockeler et al, to be published. 

[6] C. W. Bernard and M. F. L. Golterman, Phys. Rev. D 49, 486 (1994); S. R. Sharpe and N. Shoresh, Phys. Rev. D 62, 
094503 (2000). 

[7] G. Martinelli et al, Nucl. Phys. B 445, 81 (1995); M. Gockeler et al, Nucl. Phys. B 544, 699 (1999). 

[8] K. G. Chetyrkin, Phys. Lett. B 404, 161 (1997); J. A. M. Vermaseren, S. A. Larin and T. van Ritbergen, Phys. Lett. B 
405, 327 (1997); K. G. Chetyrkin and A. Retey, Nucl. Phys. B 583, 3 (2000). 

[9] A. Ali Khan et ai, Nucl. Phys B 689, 175 (2004). 
[10] C. Aubin et ai, | arXiv:hep -lat/0402030 
[11] C. R. Allton et ai, Phys. Rev. D 65, 054502 (2002). 



mf s (2GeV) 
mfS(2 GeV) 



4.7(2)(3) MeV, 
119(5)(8) MeV, 



(7) 
(8) 



